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Abstract. The paper presents mathematical models of quasicrystals with particu- 
lar attention given to cut-and-project sets. We summarize the properties of higher- 
dimensional quasicrystal models and then focus on the one-dimensional ones. For 
the description of their properties we use the methods of combinatorics on words. 
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1. Introduction 

Crystals have been admired by people since long ago. Their geometrical shape distin- 
guished them from other solids. Rigorous study of crystalline structures started in years 
1830-1850 and was crowned around 1890 by famous list of Fedorov and Schoenflies 
which contained 32 classes of crystals. Their classification was based purely on geom- 
etry and algebra. The simplest arrangement, arrangement found in natural crystals, is a 
simple repetition of a single motive. In mathematics, it is described by the lattice theory, 
in physics, the subject is studied by crystallography. Repetition of a single motive means 
periodicity. Another remarkable property, characteristic for crystals, is their rotational 
symmetry, i.e. invariance under orthogonal linear transformations. 

Important consequence of the lattice theory is that neither planar nor space (three- 
dimensional) periodic structures can reveal rotational symmetry of order 5 or higher than 
6, see Il27l . 

The discovery made by Max von Laue in 1912 enabled the study of the atomic 
structure of crystals via X-ray diffraction patterns and, in fact, justified the theoretical 
work developed by crystallography before. In case of periodic structures, the type of 
rotational symmetry of the crystal corresponds to the type of rotational symmetry of the 
diffraction diagram. 

The discovery that rapidly solidified aluminium-manganese alloys has a three- 
dimensional icosahedral symmetry, made by Shechtmann et al. L28J in 1982, was there- 
fore an enormous surprise for the crystallographic society. The diffraction diagram of 
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this alloy revealed five-fold rotational symmetry. Materials with this and other crystal- 
lographically forbidden symmetries were later fabricated also by other laboratories by 
different technologies. They started to be called quasicrystals. 

Schechtman's discovery shows that periodicity is not synonymous with long-range 
order. The definition of long-range order is however not clear. By this term crystallo- 
graphers usually understand ordering of atoms in the material necessary to produce a 
diffraction pattern with sharp bright spots. This is also used in the general definition of 
crystal adopted by Crystallographic Union at its meeting in 1992. 

The only clear requirement agreed upon by all scientists is that the set modeling qua- 
sicrystal, i.e. positions of atoms in the material, should be a Delone set. Roughly speak- 
ing, this property says that atoms in the quasicrystal should be 'uniformly' distributed in 
the space occupied by the material. Formally, the set S C K'' is called Delone if 

(i) (uniform discreteness): there exists ri > such that each ball of radius ri contains 

at most one element of S; 

(ii) (relative density): there exists r2 > such that each ball of radius r2 contains at least 

one element of S. 

The requirement of the Delone property is however not sufficient, for, also positions 
of atoms in an amorphous matter form (a section of) a Delone set. Therefore Delone 
sets modeling quasicrystals must satisfy other properties. According to the type of these 
additional requirements there exist several approaches to quasicrystal definitions |2i] 
1221 : The concept of Bohr (Besicovich) almost periodic set, is based on Fourier analysis. 
The second concept of Patterson sets is based on a mathematical analogue of X-ray 
diffraction, developed by Hof. The third concept, developed by Yves Meyer, is based on 
restriction on the set E — S of interatomic distances. It is elegant and of purely geometric 
nature: A Meyer set S e M'* is a Delone set having the property that there exists a finite 
set F such that 

E-S c E + F. 

In 1211 . Lagarias has proven that a Meyer set can equivalently be defined as a Delone set 
E such that E — E is also Delone. 

There exists a general family of sets E that are known to have quasicrystalline prop- 
erties: the so-called cut and project sets, here abreviated to C&P sets. Various subclasses 
of these sets appear to satisfy all three above quoted definitions of quasicrystals. 

The paper is organized as follows: The construction of quasicrystal models by cut 
and projection is introduced in Section|2]and illustrated on an example of cut-and-project 
set with 5-fold symmetry in Section 13 The remaining part of the paper focuses on the 
properties of one-dimensional cut-and-project sets. Section |3 provides their definition 
Sections|5J|5]and0show their diverse properties, both geometric and combinatorial. 

2. Cut-and-project Sets 

The construction of a cut-and-project set (C&P sets) starts with a choice of a full rank 
lattice: let xi , 0:2 , . . . , a;^ e M'' be linearly independent vectors over M, the set 
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L = {aixi + 02X2 + . . . + adXd \ ai, . . . , e Z} . 

is called a lattice. It is obvious that a lattice is a Delone set. Mathematical model for ideal 
crystal (or perfect crystal) in is a set A, which is formed by a finite number of shifted 
copies of a single lattice L. Formally, A is a perfect crystal if A = L + S*, where S" is a 
finite set of translations. 

Since lattices satisfy L — L C L and a perfect crystal satisfies A — A C A — 5, they 
are both Meyer sets. The Meyer concept of quasicrystals thus generalizes the classical 
definition of crystals. Perfect crystal is however a periodic set, i.e. A + a; C A for any 
X G L, and therefore it is not a suitable model for quasicrystalline materials, which reveal 
rotational symmetries incompatible with periodicity. We shall now describe a large class 
of Meyer sets which are not invariant under translation. 

Let L be a full rank lattice in R'' and let R'' be written as a direct sum Vi © V2 of 
two subspaces. One of the subspaces, say Vi, plays the role of the space onto which the 
lattice L is projected, we call Vi the physical space, the other subspace, V2, determines 
the direction of the projection map. V2 is called the inner space. Let us denote by vri the 
projection map on Vi along V2, and analogically for tt2- 

We further require that the full rank lattice is in general position, it means that tti is 
one-to-one when restricted to L, and that the image of the lattice L under the projection 
7r2 is a set dense in ¥2 - The situation can be diagrammed as follows; 



U 
L 

For the definition of CSzP sets we need also a bounded set E V2, called accep- 
tance window, which realizes the "cut". The ChP set is then defined as 

E(f2) {t^i{x) \ x E L and 7r2(x) G f2} . 

A cut-and-project set with acceptance window Vl is formed by lattice points pro- 
jected on Vi, but only those whose projection on V2 belongs to 51, i.e. by projections of 
lattice points found in the cartesian product Vi x VL. 

Figure [2 shows the construction of a CSzP set with one-dimensional physical and 
one-dimensional inner space. The acceptance window here is an interval in V2. And the 
cylinder Vi x is a strip in the plane. 

Let us list several important properties of C&P sets: 

• + t ^ for any t € Vi, i.e. T,{VL) is not an ideal crystal. 

• If the interior 51° of the acceptance window is not empty, then £(57) is a Meyer 
set^ 

• If S is a Meyer set, then there exists a C&P set I](51) with f2° 7^ and a finite 
set F such that E C E(51) + F. 

First two properties can be derived directly from the definition of a C&P set; the 
third one has been shown in 1 22 1 . 

The aim of physicists is to find a mechanism which forces atoms in the quasicrys- 
talline matter to occupy given positions. All physical approaches for describing crystals 
are based on minimum energy argument. If one wants at least to have a chance to find a 
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Figure 1. Construction of a one-dimensional cut-and-project set. 



physical explanation why a given Delone set is a suitable model for a quasicrystal then 
the number of various neighborhoods of atoms (of points) in the Delone set must be fi- 
nite. This requirement is formalized in the notion of finite local complexity: We say that 
a Delone set S is of finite local complexity if for any fixed radius r all balls of this radius 
r contain only finitely many different configurations of points of S up to translation. 

It follows from their definition, that Meyer sets have finite local complexity. There- 
fore the condition f2° 7^ ensures that a C&P set I](f2) has finite local complexity. 

Another physically reasonable requirement on the model of quasicrystal is that every 
configuration of points are found in the modeling set infinitely many times. This property 
may be for example ensured by the requirement that the boundary of the acceptance 
window has an empty intersection with the image of the lattice by the projection, i.e. 

ar2n7r2(L) = 0. 



3. Cut-and-project Sets with Rotational Symmetry 

Recalling the motivation for introducing the notion of quasicrystals, one should ask about 
conditions ensuring existence of a crystallographically forbidden symmetry. For this, 
conditions on the acceptance window alone are not sufficient. The construction of a two- 
dimensional C&P set with 5-fold symmetry has been described by Moody and Patera 
in |25|. In [6J one can find a more general construction of C&P sets with rotational 
symmetry of order 2n + 1, for n G N, n > 2. 

Consider the lattice i C K** to be generated by unit vectors ai, . . . , q;4 whose mu- 
tual position is given by the following diagram. 



ai 012 0^3 Q!4 
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In the diagram the vectors connected by an edge make an angle 7r/3, otherwise are or- 
thogonal. Such vectors are root vectors of the group which is in physics denoted by A4. 
It can be verified that the lattice generated by ai, . . . , q;4 is invariant under 5-fold rota- 
tional symmetry. Let us mention that dimension 4 is the smallest which allows a lattice 
with such a rotational symmetry. 

The physical space Vi and the inner space V2 in our example have both dimension 
2, thus they are spanned by two vectors, say v, u and u*, u* respectively. We can choose 
them as unit vectors, such that u and v form an angle ^ and u* a v* form an angle 
The following scheme shows the definition of the projection tti, which is uniquely given 
if specified on the four basis vectors ai , . . . , 0:4. 

a2 Ct4 TV V 

ai u TU 

The irrational number r, usually called the golden ratio, is the greater root of the 
quadratic equation = x + 1. Recall that regular pentagon of side-length 1 has the 
diagonal of length r, whence the correspondence of the golden ratio with the construction 
of a point set having 5 -fold rotational symmetry. 

The projection 1:2 is defined analogically, substituting vectors u, v in the diagram by 
v*,u*, and the scalar factor t by t' = ^(1 — V5), which is the other root of the quadratic 
equation x'^ = a; + 1. 

With this choice of projections tti and 7T2, a point of a lattice given by four integer 
coordinates (a, b, c, d) is projected as 

TTi (a, 6, c, d) = (a + Tb)v + (c + Td)u , 
7r2(a, 6, c, d) — {a + T'b)v* + (c + Td')u* , 

and the C&P set has the form 

= {(a + Tb)v + (c + Td)u I a,b,c,de Z, {a + T'b)v* + (c + Td')u* € f)} . 

To complete the definition of the C&P set S(ri) we have to provide the acceptance 
window n. Its choice strongly influences geometrical properties of I](r2). In 1121 it is 
shown that with the above cut-and-project scheme the set has 10-fold rotational 

symmetry if and only if the 10-fold symmetry is displayed by the acceptance window il. 
Figure|2]shows a cut-and-project set I](f2) where is a disc centered at the origin. 

If one studies the inter-atomic interactions, it is impossible to consider contribution 
of all atoms in the matter; one must limit the consideration to 'neighbours' of a given 
atom. Thus it is necessary to define the notion of neighbourhood in a general point set, 
which has not a lattice structure. A natural definition of neighbours is given in the notion 
of a Voronoi cell. 

Consider a Delone set E C M"* and choose a point x in it. The Voronoi cell of x is 
the set 



where r = . 
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V{x) = {y e M'* I \\x-v\\ < \\z-y\\ for all z e E} . 

The Voronoi cell of the point x is thus formed by such part of the space, which is 
closer to x than to any other point of the set S. Since S is a Delone set, the Voronoi cells 
of all points are well defined convex polytopes in K'', filling this space without thick 
overlaps. Therefore they form a perfect tiling of the space. 

The notion of Voronoi cells allows a natural definition of neighbourhood of points in 
a Delone set S C M''. Two points may be claimed neighbours, if their Voronoi polytopes 
share a face of dimension d — 1. The Voronoi tiling of the cut-and-project set from 
Figure|2]is shown on Figure|3l 

In the Voronoi tiling of Figure |3l one finds only 6 basic types of tiles (Voronoi 
polygones). They are all found together with their copies rotated by angles j^j, j = 
0,1,. ..,9. 

Let us mention that the Voronoi tiling shown at Figurel^is aperiodic, since the C&P 
set I](f2) is aperiodic. For certain classes of acceptance windows with 10-fold rotational 
symmetry, the collections of appearing Voronoi tiles are described in the series of arti- 
cles |23|. 

The geometry of the Voronoi tilings generated by cut-and-project sets is, except 
several special cases, unknown. The only known fact is that the number of types of tiles 
is for every cut-and-project set finite, which follows from the finite local complexity of 
cut-and-project sets. The situation in dimenions d > 2 is quite complicated. Therefore 
we focus in the remaining part of the paper on one-dimensional cut-and-project sets. 



4. One-dimensional CkP Sets and C&P Words 



Let us describe in detail the construction of a one-dimensional Ck,P set, as illustrated 
on Figure^ Consider the lattice L = 1? and two distinct straight lines Vi : y = ex and 
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Figure 3. Voronoi tiling of the cut-and-project set shown in Figurel2l 

^ i ^ 6 O 

Figure 4. The tiles appearing in the Voronoi tihng of Figurel3l 

^2 : 2/ = V^, e ^ rj.lf we choose vectors xi — and X2 — ^37(1, ?7) then 

for any point of the lattice we have 

{p, q) = {q - pf])xi + {q - pe)x2^ . 

Let us recall that the construction by cut and projection requires that the projection tti 
restricted to the lattice L is one-to-one, and that the set Ti2{L) is dense in V2. 

If 77 and e are irrational numbers, then these conditions are satisfied. The projection 
of the lattice L = on the straight lines Vi and V2 are written using additive abelian 
groups 

:={a + 677 I a,6e Z}, 
Z[£] :={a + 6e I a,6e Z}. 

These groups are obviously isomorphic; the isomorphism ★ : Z[r;] Z[e] is given by 
the prescription 

X = a + bri i—t X* = a + be . 



The cut-and-project scheme can then be illustrated by the following diagram. 



g 
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Every projected point 7ri(p, (7) lies in the set lJ\r^x\. The notation will be simplified 
by omitting the constant vector x\. In a similar way, we omit the vector in writing 
the projection ■n2{p^ <?)■ With this convention, one can define a one-dimensional cut-and- 
project set as follows. 

Definition 1. Let e, be distinct irrational numbers, and let c M be a bounded inter- 
val. Then the set 

{a + 6?7 I a, 6 e Z, a + e f7} = {x G 7L\i]\ | e f7} 
is called a one-dimensional cut-and-project set with parameters e, 77 and il. 

From the properties of general cut-and-project sets, in particular from their finite 
local complexity, we derive that the distances between adjacent points of I]j_,,(ri) are 
finitely many. The following theorem 1 19 1 limits the number of distances to three. 

Theorem 2. For every Eg^^(r2) there exist positive numbers Ai, A2 S ^[77] such that 
the distances between adjacent points in ,;(ri) take values in {Ai, A2, Ai + A2}. 
The numbers Ai, A2 depend only on the parameters e, rj and on the length of the 
interval $1. 

As the theorem states, the distances Ai , A2 and Ai + A2 depend only on the length 
of the acceptance window Q. and not on its position in R or on the fact whether is open 
or closed interval. These properties can, however, influence the repetitivity of the set 
Eg ,,(^1). More precisely, they can cause that the largest or the smallest distance appears 
in (il) at one place only. 

From the proof of Theorem|2lit follows that if VI semi-closed, then E^ (il) is repet- 
itive. Therefore in the sequel we consider without loss of generality the interval Q. of the 
form il — [c, c + ^) . Nevertheless, let us mention that even if Ee_^ ( [c, c + ) is repetitive, 
the distances between adjacent points may take only two values, both of them appearing 
infinitely many times. 

The knowledge of values Ai, A2 allows one to easily determine the neighbours of 
arbitrary point of the set Ee_,,([c, c + €)) and by that, generate progressively the entire 
set. Denote {xn)nez an increasing sequence such that {x„ | n G Z} = E£_^([c, c + £)). 
Then for the images of points of the set E^ ,j([c, c + £)) under the mapping ★ one has 

if < e [c,c + £-A*^) =: Qa, 
if <e[c + ^-A^c-A^)=:f]B, (1) 
if < e [c- A5,c + £) =: Qc- 

The mapping, which to associates the element x*^_^^i is a piecewise linear bijection 
/ : [c,c + i) 1^ [c,c + £). Its action is illustrated on Figure |5l Such mapping is in 
the theory of dynamical systems called 3-interval exchange transformation. In case that 



X* + A{ + AJ 
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A{ — A2 — \^\ = £, the interval denoted by fi^ is empty. This is the situation when the 
distances between adjacent points in Eg ([c, c + £)) take only two values. The mapping 
/ is then an exchange of two intervals. 




Figure 5. The diagram illustrating the prescription 0- The function / is a three-interval exchange transfor- 
mation. 

In order to record some finite segment of the set I]£,,j([c, c + £)), one can write 
down to a list the individual elements of this set. However, this is not the most efficient 
way, since every point x G I]£^^([c, c + £)) is of the form x ~ a + brj E Z[ri] and 
with growing size of the considered segment, the integer coordinates a, b, of the points 
needed for recording x = a + bij grow considerably. Much simpler is to find a point 
xq E I]£,,,([c, c + £)) in the segment and record the sequence of distances between 
consecutive points on the left and on the right of xq. 

With this in mind, the entire set !]£,,,( [c, c + £)) can be coded by a bidirectional 
infinite word (it„)„gz in the alphabet {A, B, C} given as 



A 


if 




Xn 


= Ai, 


B 


if 


Xn+l 


Xfi 


= Al + A2 


C 


if 


Xn+l 


Xn 


= A2. 



Such an infinite word is denoted by 

Example 1 (Mechanical words). Let us choose irrational e e (— 1, 0) and irrational 77 > 
0. We shall consider one-dimensional C&P set with acceptance window VI — (/3 — 1, 
for some /3 G K. For simplicity of notation we shall put a = —e E (0, 1). From the 
definition of C&P sets it follows that 

a + bri eT,^^^{Q.) ^ a + beEfl <^ (3-l<a-ba<P ^ a^[ba + P\ 
and therefore the CSzP set is of the form 

S„„,^(^ -l,f3]^{lba + l3\+bf]\bEZ}. 

Since a, 77 > 0, the sequence x„ := [na + P\ + nrj is strictly increasing and thus the 
distances between adjacent points of the C&P set S_a ,,(/3 — 1, /3] are of the form 

Xn+i-Xn = V+ [{n + l)a + l3\ - [na + P\ = 1'^ J^l ' 
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This ChP set can therefore be coded by an infinite word (u„)„gz in a binary alphabet. 
Usually one chooses the alphabet {0, 1}, so that the n-th letter of the infinite word can 
be expressed as 

Un= [{n+l)a + I3\- \iia + I3\ . (3) 

Such infinite words {un)ne7. were introduced already in 1261 and since then, extensively 
studied, under the name mechanical words. The parameter a G (0, 1) is called the slope 
and the parameter (3 the intercept of the mechanical word. 

5. Equivalence of One-dimensional C&P Sets 

In the previous section we have defined for a triplet of parameters e, 77 and the set 
Eg and we have associated to it the infinite word Ue.^(r2). Natural question is how 
these objects differ for different triplets of parameters. Recall that our construction is 
based on the projection of the lattice 1?. The group of linear transformations of this 
lattice onto itself is known to have three generators: 




Directly from the definition of C&P sets it follows that the action of these three genera- 
tors on the lattice provides the identities 

=77El 

Another identity for C&P sets is obtained using invariance of the lattice 7? under 
translations. 

a + b-q + Eg ,,(51) = ^^^^(ri + a + be) , for any a,b eZ . 

The mentioned transformations were used in fT9\ for the proof of the following 
theorem. 

Theorem 3. For every irrational numbers e^rj, e =/= rj and every bounded interval f2, 
there exist e, f) and an interval 17, satisfying 

(P) £€(-1,0), 7^>0, max(l + £,-e) < li^l < 1 

such that 



Eg^^(il) — sEe_^(il) , for some s e M. 
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Multiplying the set ,^(17) by a scalar s can be understood as a choice of a different 
scale. From the physical point of view the sets are therefore de facto the same. In geom- 
etry such sets are said to be similar In the study of one-dimensional C&P sets, one can 
therefore limit the considerations only to parameters satisfying the condition [P). One 
may ask whether the family of parameters can be restricted even more. More precisely, 
one asks whether for different triples of parameters satisfying (P) the corresponding 
C&P sets are essentially different. The answer to such question is almost always affir- 
mative, except certain, in some sense awkward, cases. A detailed analysis can be found 
in Gil- 

Theorem|3lconcerned geometrical similarity of C&P sets. If interested only in the 
corresponding infinite words Ue_,,(r2), we can restrict the consideration even more. This 
is a consequence of the following two assertions. 

Claim 1. If the parameters £,rji^Q, and e,r]2,^ satisfy (P) then the infinite word 
Ue^nii^) coincides with u^j-i^iVl). 

Consequently, we can choose the slope of the straight line V2 in the cut-and-project 
scheme to be 77 — — i, where e is the slope of the straight line Vi. The straight lines Vi 
and V2 can therefore be chosen without loss of generahty mutually orthogonal. 

Claim 2. If the parameters e, 77, VL satisfy (P) then the infinite word Ue,r]{^) coincides 
with 17) up to permutation of assignment of letters. 

This statement implies that for the study of combinatorial properties of infinite words 
associated with C&P sets, one can limit the choice of the parameter e to the range 

(-io). 

6. Factor and Palindromic Complexity of C&P Words 

For the description of combinatorial properties of infinite words associated to one- 
dimensional C&P sets one uses the terminology and methods of language theory. Con- 
sider a finite alphabet A and a bidirectional infinite word u = {un)nez, 

U — ■ ■ ■ U^2U-iUqUiU2 ■ ■ ■ ■ 

The set of factors of u of the length n is denoted 

Cn = {UiUi+i ■ ■ ■ Ui+n-l I i e Z} . 

The set of all factors of the word u is the language of u. 

If any factor occurs in u infinitely many times, then the infinite word u is called recurrent. 
If moreover for every factor the gaps between its individual occurrences in u are bounded, 
then u is called uniformly recurrent. 

The factor complexity of the infinite word m is a mapping C : N ^ N such that 
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C[n) = • • • Ui+n-l I « G = #£„ . 

The complexity of an infinite word is a measure of ordering in it: for periodic words it is 
constant, for random words it is equal to 

Since every factor UiU^+i • • • Ui+n-i of the infinite word u = (u„)„gz has at least 
one extension UiUi+i ■ ■ -Ui+n, it is clear that C{n) is a non-decreasing function. If 
C{n) — C(n + 1) for some n, then every factor of the length n has a unique extension 
and therefore the infinite word u is periodic. The complexity of aperiodic words is nec- 
essarily strictly increasing function, which implies C{n) > n + 1 for all n. It is known 
that mechanical words (defined by (|3}) with irrational slope are aperiodic words with 
minimal complexity, i.e. C{n) = n + 1. Such words are called bidirectional sturmian 
words. A survey of which functions can express the complexity of some infinite word 
can be found in 1 16,1. 

The following theorem has been proven in II19I for the infinite words obtained by the 
cut-and-project construction. 

Theorem 4. Let Ue^,,(ri) be a CSzP infinite word with £7 = [c, c + t). 

• If t. ^ then for any n g N we have C{n) = 2n + 1. 
» If i € Z,[e], then for any n G N we have C{n) < n + const. 

Moreover, the infinite word ,,(17) is uniformly recurrent. 

From the theorem it follows that the complexity of the infinite word (fl) depends 
only on the length ^ = of the acceptance window and not on its position. This is the 
consequence of the fact that language of depends only on \ Vl\. 

If the parameters e, -q and VL = [c, c+£) satisfy the condition (P), then the complexity 
of the infinite word is minimal (i.e. C{n) = n + 1) if and only if ^ = 1. Thus the 

infinite words [c, c + 1)) are sturmian words. Nevertheless, the sturmian structure 
can be found also in words Ue.,,([c, c + 1^)) with I e Z[e]. For, one can prove the 
following proposition [ISJ . 

Proposition 5. If £ G Z[e] then there exists a sturmian word v = ■ ■ ■ t)_2U-ii'oi'ii'2 • • • 
in {0, 1}^ and finite words Wq, Wi over the alphabet {A, _B, C} such that 

U,^r,{^) = ■■■ W,^,W,_,W^,W^,W^, ■■■ . 

The proposition in fact states that the infinite word ,,(17) can be obtained by con- 
catenation of words Wo, W^i in the order of O's and I's in the sturmian word v. Let us 
mention that Cassaigne [11 1 has shown a similar statement for arbitrary one-directional 
infinite words with complexity n + const. He calls such words quasisturmian. 

A reasonable model of quasicrystalline material cannot distinguish between the or- 
dering of neighbours on the right and on the left of a chosen atom. In terms of the in- 
finite word, which codes the one-dimensional model of quasicrystal, it means that the 
language C must contain, together with every factor w = wiW2 . ■ ■ Wn also its mirror 
image w — WnWn-i . . . wi. The language of ChP words satisfies such requirement. 

A factor w, which satisfies w = w, is called a palindrome, just as it is in natural 
languages. The study of palindromes in infinite words has a great importance for de- 
scribing the spectral properties of one-dimensional Schrodinger operator, which is as- 
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sociated to (un)ngz in the following way: To every letter of the alphabet a ^ A one 
associates the potential V{a) in such a way that the mapping y : ^ i-^ M is injective. 
The one-dimensional Schrodinger operator H is then defined as 

{H4)){n) = (f){n + 1) + (f){n - 1) + V{un)(t>{n) , (t> E . 

The spectral properties of H influence the conductivity properties of the given struc- 
ture. Roughly speaking, if the spectrum is absolutely continuous, then the structure be- 
haves like a conductor, while in the case of pure point spectrum, it behaves like an insu- 
lator. In | 20J one shows the connection between the spectrum of H and the existence of 
infinitely many palindromes in the word (w„)„gz. 

The function that counts the number of palindromes of a given length in the lan- 
guage C of an infinite word u is called the palindromic complexity of u. Formally, the 
paUndromic complexity of u is a mapping : N N defined by 

V{n) := G £„ | w = w} . 

Upper estimates on the number 7^(71) of palindromes of length n in an infinite word 
u can be obtained using the factor complexity C{n) of u. In |3| the authors prove a 
result which puts in relation between the factor complexity C{n) and the palindromic 
complexity V^n). For a non-periodic infinite word u it holds that 

V{n)<-C(n+\^\] . (4) 

Combination of the above estimate with the knowledge of the factor complexity we 
obtain for CkP infinite words that Tin) < 48. 

Infinite words constructed by cut and projection are uniformly recurrent. For such 
words, the upper estimate of the palindromic complexity can be improved, using the 
observation that uniformly recurrent words have either — for sufficiently large 

n, or the language C of the infinite word is invariant under the mirror image, see [5|. If 
C contains with every factor w its mirror image w, then 

V{n) + V{n + 1) < 3AC(n) i{c{n + 1) - C(n)) . (5) 

This estimate of the palindromic complexity is better than that of in case that the fac- 
tor complexity C{n) is subpolynomial. In particular, for C&P infinite words we obtain 
'P{n) < 6. In |il4| one can find the exact value of the palindromic complexity for infi- 
nite words coding three-interval exchange transformation. Since this is the case of C&P 
words, we have the following theorem. 

Theorem 6. Let Ue^r/i^) be a C&P infinite word with VL — [c, c + ^) and let e, rj, fl 
satisfy the conditions (P). Then 

{1 for n even, 
2 for n odd and £ ~ 1, 
3 for n odd and t <\. 
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7. Substitution Invariance of CkP Words 

To generate the set using the definition resumes in deciding for every point of 

the form a + br], whether a + he belongs to the interval Q, or not. This is done by verifying 
certain inequalities between irrational numbers. If we use a computer working with finite 
precision arithmetics, the rounding errors take place and in fact, the computer generates 
a periodic set, instead of aperiodic ,,(17). The following example gives a hint to much 
more efficient and in the same time exact generation of a C&P set. 

Consider the most popular one-dimensional cut-and-project set, namely the Fi- 
bonacci chain. It is a C&P set with parameters rj = t, e = t' and — [0, 1). (Recall 
that the golden ratio r = ^^/^ and r' = are the roots of the equation — x+1.) 

Since = t + 1, the set of all integer combinations of 1 and r is the same as the set of 
all integer combinations of and r, formally 

tZ[t] = Z[r] . 

Moreover, Z[r] is closed under multiphcation, i.e. Z[t] is a ring. Since r + r' = 1, we 
have also Z[t] = Z[r'], and the mapping -k which maps a + br a + br' is in fact an 
automorphism on the ring Z[t] . Note that Z[t] is the ring of integers in the field Q[t] and 
* is the restriction of the Galois automorphism of this field. 

Using the mentioned properties one can derive directly from the definition of C&P 
sets that 

which is valid for every acceptance window ft. In the case J7 = [0, 1) we moreover have 
(T')^n C n. Therefore 

i.e. Yir'.ri^) is selfsimilar with the scaling factor r^, as illustrated on Figure|S] 

Example [2 namely equation (|3 implies that Yir'.ri^) has two types of distances 

between adjacent points, namely Ai = and A2 = t. In Figure |6lthe distance Ai is 

coded by the letter A and the distance A2 by the letter B. 

For our considerations it is important that every distance A scaled by the factor is 

filled by two distances A followed by B. Similarly, every scaled distance B is filled by 

A followed by B. This property can be proven from the definition of E^' ,-(17). 

As a consequence, the Fibonacci chain can be generated by taking an initial segment 

of the set, scaling it by and filling the gaps by new points in the above described way, 

symbolically written as the rule 

A 1-^ AAB and B AB. 

Repeating this, one obtains step by step the entire C&P set. Since the origin as an 
element of I]^/^^(ri) has its left neighbour in the distance A2 and the right neighbour in 
the distance Ai, we can generate the Fibonacci chain symbolically as 



B\A ^ AB\AAB ^ AABAB\AABAABAB ^ 
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Figure 6. Selfsimilarity and substitution invariance of tlie Fibonacci word. 

A natural question arises, whether such efficient and exact generation is possible also 
for other one-dimensional cut-and-project sets, respectively their infinite words. Let us 
introduce several notions which allow us to formalize this question. 

A mapping ip on the set of finite words over the alphabet ^ is a morphism, if the ip- 
image of a concatenation of two words is concatenation of the 9?-images of the individual 
words, i.e. (p{vw) = ip{v)(f{w) for every pair of words v, w over the alphabet A. For 
the determination of a morphism, it suffices to specify the (^-images of the letters of the 
alphabet. The action of a morphism can be naturally extended to infinite words, 

(p{u) = </9(. . . U-2U-i\uoUiU2Us ...):=... ip{U-2)ip{U-i)\ip{ui)(p{U2)(p{U3) . . . 

An infinite word u invariant under the action of the morphism, i.e. which satisfies ip{u) = 
u, is called a fixed point of (p. In this terminology, one can say that the Fibonacci chain (or 
the infinite word coding it) is a fixed point of the morphism over a two-letter alphabet 
{A,B}, which is determined by the images of letters, ip{A) = AAB, ip{B) = AB. 

The identity map, which maps every letter of the alphabet on itself, is also a mor- 
phism and arbitrary infinite word is its fixed point. However, one cannot use the identity 
map for generation of infinite words. Therefore we must put additional requirements on 
the morphism (p. 

The morphism ip over the alphabet A is called a substitution, if for every letter a £ A 
the length of the associated word ip{a) is at least 1, and if there exist letters ao, bo & A 
such that the words (p{ao) and (p{bo) have length at least 2, the word (p{ao) starts with 
the letter ao, and the word </?(6o) ends with the letter bo- 

A morphism, which is in the same time a substitution, necessarily has a fixed point 
u, which can be generated by repeated appUcation of the morphism on the pair of letters 
6o|cio- Formally, 

ip{u) = u= lim .^"(6o)|¥'"(ao) • 

n — *oo 

To every substitution (p over the alphabet A = {ai, a2,. ■ . ,ak} one associates a 
k X k square matrix M (the so-called incidence matrix of the substitution). The element 



16 



E. Pelantovd et al. / Quasicrystals 



Mij is given as the number of letters aj in the word (p{ai). The incidence matrix of the 
substitution generating the Fibonacci word is M — (il) ■ 

The incidence matrix of a substitution is by definition a non-negative matrix for 
which Perron-Frobenius theorem holds 1 17.1 . A substitution ip is said primitive, if some 
power of its incidence matrix is positive. In this case the spectral radius of the matrix 
is an eigenvalue with multiplicity 1, the corresponding eigenvector (the so-called Perron 
eigenvector of the matrix) being also positive. 

Although the incidence matrix M does not allow to reconstruct the substitution (p, 
many properties of the fixed points of ip can be derived from it. Let us mention some of 
them. 

• If the infinite word u is invariant under a substitution, then there exists a constant 
K such that for the complexity function of the word u we have 



C{n) < K-r? for all n e N. 



• If the infinite word u is invariant under a primitive substitution, then there exist 
constants K\ and K2 such that for the factor complexity and palindromic com- 
plexity of the word u we have 



On) < Kin and V{n) < K2 for all n e N . 



• An infinite word which is invariant under a primitive substitution is uniformly 
recurrent. 

• If the infinite word u is invariant under a primitive substitution Lp over an alphabet 
A = {ai, 02, ... , Cfc}, then every letter has well defined density in u, i.e. the 
limit 

number of letters in the word u_„ ... u_i|uoUi ... Un_i 
p{a,) := hm — — 

rwoo Zn + I 

exists. Let {xi, X2, ■ ■ ■ Xk) be the Perron eigenvector of the matrix (transpose 
of the incidence matrix M of the substitution p). Then the density p(ai) is equal 
to 



Xi + X2 + . ■ ■ + Xk 



The question of describing all CSzP infinite words invariant under a substitution is 
stiU unsolved. A complete solution is known only for CSzP words over a binary alpha- 
bet, which can be, without loss of generality, represented by mechanical words ^ with 
irrational slope a G (0, 1) and intercept (3 G [0, 1). 

The substitution invariance of mechanical words has first been solved in 1131 for the 
so-called homogeneous mechanical words, i.e. such that (3 — Q. 

Theorem?. The homogeneous mechanicalword with slope a € (0,1) is invariant under 
a substitution if and only if a is a quadratic irrational number whose conjugate a' does 
not belong to (0, 1). 
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A quadratic irrational number a E (0, 1) whose conjugate a' ^ (0, 1) is called 
Sturm number. Let us mention that in the paper 1131 the Sturm number is defined using 
its continued fraction expansion. The simple algebraic characterization was given in |2|. 

The substitution invariance for general (inhomogeneous) mechanical words is solved 
in|4|and(29|. 

Theorem 8. Let a be an irrational number, a e (0, 1), and let f3 e [0, 1). The mechani- 
cal word with slope a and intercept (3 is invariant under a substitution if and only if the 
following three conditions are satisfied: 

(i) a is a Sturm number, 

(ii) 13 e Q[a], 

(iii) a' < (3' < 1 — a' or 1 — a' < (3' < a' , where [3' is the image of [3 under the 
Galois automorphism of the field Q[q!]. 

Unlike the case of binary C&P words, the question of substitution invariance for 
ternary C&P words has been sofar solved only partially. The following result is the 
consequence of I1I9I . 

Theorem 9. Let = [c, d) be a bounded interval. If the infinite word We^,,(r2) is in- 
variant under a primitive substitution, then e is a quadratic irrational number and the 
boundary points c, d of the interval belong to the quadratic field Q(£). 

All C&P words satisfying the properties of the theorem have a weaker property than 
substitution invariance, the so-called substitutivity 1 15|. It allows one to generate even 
those infinite words which are not fixed points of a morphism. 

8. Conclusions 

In the theory of mathematical quasicrystals, best known are the properties of the one- 
dimensional models, be it the geometric or the combinatorial aspects of these structures. 
However, this information can be used also in the study of higher dimensional models, 
since the one-dimensional ones are embedded in them. In fact, every straight line con- 
taining at least two points of a higher-dimensional cut-and-project set, contains infinitely 
many of them, and they ordering is a one-dimensional cut-and-project sequence. 

Nevertheless, the notions of combinatorics on words, as they were presented here, 
are being generalized also to higher dimensional structures; for example, one speaks 
about complexity and substitution invariance of two-dimensional infinite words, even 
two-dimensional sturmian words are well defined [7 8 1. 

Except cut-and-project sets, there are other aperiodic structures which can serve 
for quasicrystal models; they are based on non-standard numeration systems [101. The 
set of numbers with integer /^-expansions share many properties required from one- 
dimensional quasicrystal models, in particular, they are Meyer sets, are self-similar, and 
the corresponding infinite words are substitution invariant. 
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